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Abstract 

Axion fields provide the most elegant solution to the strong CP problem. In 
string compactifications it is difficult to obtain an axion whose decay constant is 
consistent with current cosmological bounds. We examine this question in theories 
with accidental U(l) symmetries that appear as low energy remnants of discrete 
symmetries. We refer to the axion-like particles from the spontaneous breakdown of 
these symmetries as accions. In such systems, the accion decay constant depends on 
the vacuum configuration and can be lowered to fit the bounds. Furthermore, we 
find that such accions with consistent decay constant can be embedded in special 
vacua of Zg-II orbifold models with realistic features. 



1 Introduction 



The most elegant and appealing solution to the strong CP-problem is based on the conjec- 
ture of an axion field [1,2]. It requires the existence of an anomalous global Peccei-Quinn 
symmetry [3] U(1) P q and its spontaneous breakdown at a scale F a (where F a denotes the 
"axion decay constant"). Constraints (mostly) from astrophysics and cosmology require 
F a to be in the axion window 

10 9 GeV <F a < 10 12 GeV (1) 

for the so-called "invisible" axion [4,5]. The axion field adjusts its vacuum expectation 
value (VEV) to cancel the ^-parameter of quantum chromodynamics (QCD) to avoid 
CP-violation due to strong interactions. 

In the presence of gravity, the existence of an exact global U(l)-symmetry might be 
problematic. This holds explicitly in the context of string theory, where potential anoma- 
lous U(l)-symmetries appear as gauge (not global) symmetries [6]. Does this exclude the 
axion solution to the strong CP-problem in the framework of string theory? 

In the present work, we focus on a set of realistic string extensions of the minimal 
supersymmetric standard model (MSSM) in the framework of the Es x Es heterotic 
string. We address the following two questions: 

• how can we possibly obtain a global U(1) P q? 

• can we disconnect the axion scale F a from the string scale M s (which is typically 
much larger, e.g. 10 16 - 10 18 GeV)? 

The analysis of the models reveals (among the various local U(l)-symmetries) the ap- 
pearance of a multitude of discrete symmetries of various types [7,8]. In such a situation 
we might find accidental U(l)-symmetries [9] if we restrict our attention to superpoten- 
tials of a limited degree in the polynomial of the fields. Higher order non-renormalizable 
terms would break this accidental U(l) P Q-symmetry back to the exact discrete symme- 
try. In the standard model of particle physics (SM) such accidental symmetries appear 
as U(1)b,l where B(L) denote baryon(lepton) number, respectively. 

Including these accidental symmetries will allow for a multitude of axion candidates 
in the view of the first question. Apart from the strong CP-problem these fields might be 
of relevance for hidden sector gauge groups (embedded in the second Eg) or quintessential 
axions [11, 12]. 

This new paradigm of string theory, the appearance of accidental global symmetries 
from exact discrete symmetries, might find interesting applications in the creation of 

1 Similar ideas have been explored in the context of Calabi-Yau compactifications [10]. 
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hierarchical structures as explained in [13]. This might then give us the freedom to lower 
the axion scale away from the string scale and thus answering our second question. 

Within the application to the axion system, pioneering work has been done by Choi, 
Kim and Kim [14] in a model based on heterotic string constructions via the Zi2-orbifold. 
Unfortunately, however, in this model the axion decay constant of the QCD axion could 
not be decoupled from the string scale. In the present paper we analyze the models 
of the so-called heterotic Mini-Landscape [15-19] in view of potential realistic axion 
candidates. This requires a careful study of a multi- axion system including a multitude 
of global and discrete symmetries both of anomalous and non-anomalous nature. Our 
analysis of section 2 reveals the fact that the properties of the system (apart from the 
structure of the model itself) will strongly depend on the vacuum configuration. 

In section [3] we shall then discuss some representative examples of the Mini- 
Landscape and show their potential suitability for the answers of the questions raised 
above. This discussion requires extensive computer calculations and we limit ourselves 
here to the conceptual questions towards the solution of the strong CP-problem via a 
suitable axion field. The construction of a complete and fully realistic extension of the 
MSSM will be relegated to future investigations. 

We shall see that accidental symmetries that naturally appear in the framework of 
string constructions are very well suited for satisfactory axion candidates with a decay 
constant in the axion window. Section H] will summarize our results, address remaining 
problems and will give an outlook for future research. 



2 Multi- Axion Systems 

The models under consideration have, in addition to some non-anomalous U(l)g auge 
gauge symmetries, (accidental) anomalous U(l) symmetries whose breaking provides 
many axion candidates. This leads to a complicated vacuum structure that requires a 
careful investigation. Here we rely on previous work [15-19] where the supersymmetric 
vacua corresponding to D- and F-fiat solutions have been analyzed. The models typically 
have onj§ anomalous gauge symmetry U(1)a [21] with a non-trivial Fayet-Iliopoulos (FI) 
term which in turn induces non- vanishing VEVs for some of the scalar fields. These break 
some of the U(l) gauge symmetries as well as the U(1)a at a scale one or two orders of 
magnitude below the string scale M s . We now also have to analyze the fate of the possible 
accidental (anomalous) U(l) symmetries. To answer the question towards a satisfactory 
axion within the "allowed window" we have to understand not only the model but also 
the vacuum configuration, which we discuss now. 

2 More are possible in the blow-up version of the models [20] 
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2.1 Axions from a Single U(l) 



Let us start with the simplest case, a super symmetric theory with standard model (SM) 
gauge group and one additional global symmetry, U(1) P q, which is assumed to have a 
mixed U(l) P q — SU(3)c — SU(3)c anomaly. Consider the matter superfields (ipi,ipi) with 
U(1) P q charges q l . If one scalar SM singlet field, say ipo, with non-vanishing PQ charge 
attains a VEV vq, its phase field a becomes an axion and induces the following term in 
the effective Lagrangian 



a g 2 

t-pff D — tr 

F a 32vr 2 



(2) 



where g denotes the SU(3)c* coupling constant, G^ u is the color field strength, and G^ u 
its dual. The axion decay constant F a in this scenario is given by 

F. = ^ (3) 

that does not depend on the normalization of the PQ symmetry. Here A is the coefficient 
of the U(1) PQ - SU(3) C - SU(3) C anomaly 

A = ( 4 ) 

with ii being the SU(3)c quadratic index of the fermion ipi, such that ii = 1 for i^i 
transforming as a triplet. Note that the chiral PQ transformation of the axion, a — ► 
a + q° vq a (with a being the transformation parameter), reproduces as expected the 
original anomaly, 

t 
32vr 2 

It is straightforward to generalize this result to the case where various scalar fields 
obtain a VEV. If the fields tpi attain VEVs Uj, the resulting axion a is given in terms of 
the phase fields Oj by 



£ eff - £ eff + aA 4^ tr GTG^ . (5) 



1 



^ a% q % vA , (6) 



where M defines the scale at which U(1) P q is broken. The axion decay constant is then 
given by 

(T) 

Considering the PQ transformation of the phase fields a,, one recovers eq. ([5]). We would 
like to point out that in this case the axion decay constant is dominated by the largest 
VEV. This behavior might change if the low energy theory contains more than one PQ 
symmetry, as we shall discuss below. 
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2.2 Axions from Multiple U(l)s 



In a next step we consider a supersymmetric theory with two PQ symmetries, denoted 
by U(l)p x U(1)q, with respective (non- vanishing) anomalies Ap and Aq, analogously 
defined as in eq. ([3]). The scalar fields of the theory transform as 



(8) 



under these symmetries. Note that the U(l) resulting from the linear combination 
q) = -AQqp-ApqQ 



(9) 



is always anomaly-free. 

Let us suppose that both U(l)s are broken spontaneously. This occurs if, for instance, 
two fields ip\ and <p2 develop VEVs, provided that their charges (qp, qq) and (qp, qq) 
are linearly independent, i.e. 



QpQq ¥= QqQp 



(10) 



If the charges do not satisfy eq. (fTUl) . only one U(l) would be broken. As before, the 
spontaneous breaking of these U(l)s induces the term 



C 



off 



a g 
F a 32^ 2 



tr 



G^G 



«1 _«2_ \ 9 

F ai F a2 ) 32vr 2 



tr 



fyjUQ 



/(!' 



(11) 



The constants F ai are found by comparing the variation of C c s due to the chiral PQ 
transformation of the phase fields, <2j — > ai + Vi (q l P a + q l Q (3), and the one expected from 
the anomalous nature ofU(l)p x U(1)q, 



£ eS -> £ cff + (aAp + (5Aq) 
The solutions are 



32^ 2 



tr 



G^ V G 



I"' 



n 



F 



a 2 



V 2 



Qp Qq ~ Qq Q 2 P 



(13) 



Note that these constants do not depend on the basis chosen for the U(l)s. Then, the 
axion and its decay constant are respectively given by 



-ai qj v 2 + a 2 qj v\ 



Fa 



1 



+ 



Fa 2 



-1/2 



vi V2 {q P q Q - q Q q P ) 

^(Q)vi? + (q)v 2 f ' 



(14a) 
(14b) 



A few comments are in order. First, one can easily verify that, contrary to the previ- 
ous case, if two U(l)s are broken by two fields attaining VEVs, the axion decay constant 
given by eq. (|14b|) is dominated by the smaller VEV: this is the key observation in view 
of a realization of a decay constant within the allowed "axion window" . Secondly, apart 
from the axion o, clearly the theory must have a Goldstone boson. It is the combination 
orthogonal to the axion, 

ai q) vi + a 2 q 2 f v 2 

h = - =d=, (15) 

^>l) 2 + (^ 2 ) 2 

associated to the broken anomaly-free symmetry, U(l)/. If this symmetry is gauged, h is 
'eaten' via the Higgs mechanism by the gauge boson of U(l)/, making it massive. If U(l) / 
is global, h remains as a massless (uncharged) field. As a last remark, from the above 
discussion we see that in theories with multiple spontaneously broken U(l)s, one could 
readily identify the axion as the unique combination of the phase fields that is orthogonal 
to all Goldstone bosons. This observation will help us in the following discussion. 



2.3 The general case 

The above two examples are sufficient to obtain an understanding of the general situation, 
where we have more than two U(l) symmetries and where the number of SM singlets 
acquiring VEVs exceeds the number of symmetries that are broken. In order to shed 
some light on the form of the axion and its decay constant, let us examine the scenario 
in which U(l)p x U(1)q gets broken when three fields attain VEVs. One combination of 
the three phase fields en will provide the CP axion a, whereas the other two orthogonal 
combinations will yield a Goldstone boson h that does perceive chiral transformations, 
and an invariant phase field a'. 

To find the axion, we can proceed as before and determine the constants F ai . A 
more manageable approach is to identify first the fields a' and h. The axion will be then 
the combination of a, orthogonal to them. The result is a generalization of eqs. (|14p . In 
particular, assuming a large hierarchy of VEVs, for instance v\ S> v 2 ^> v%, the axion 
decay constant becomes 

F a ~ v 2 i , (16) 

If 

that is, F a is dominated by the second largest VEV (v 2 in our example) and neither by 
the largest nor by the smallest. The second largest VEV corresponds to the largest scale 
at which both symmetries, U(l)p and U(1)q, are already broken. So, it seems natural to 
expect that, for any number of SM fields attaining VEVs, the axion decay constant will 
be dominated by the scale at which all PQ symmetries get broken. In the current case, 



6 



with two PQ symmetries and a large VEV hierarchy, the scale of F a will be given by the 
second largest VEV. 

We can now state the result for the general case of M U(l) symmetries and N > M 
SM singlet fields that obtain a VEV. By Goldstone's theorem, we have generically M 
"Goldstone bosons" from the M spontaneously broken symmetries. This means that we 
have N — M invariant combinations of phase fields, which are orthogonal to these Gold- 
stone bosons. Suppose that one of the Goldstone bosons lies in the anomalous direction 
and becomes an axion, while the other M — 1 ones remain massless, possibly eaten by 
gauge bosons. From the previous discussion, if the U(l)s are broken at hierarchically 
different scales, then the axion decay constant F a is dominated by the Mth largest VEV. 

2.4 Embedding in the heterotic string theory 

In the heterotic string, there are various sources for axions. The spontaneous breakdown 
of the pseudo-anomalous U(1)a gauge symmetry together with the Green-Schwarz mech- 
anism produce the so-called model-independent axion, whose decay constant is fixed by 
the Planck scale (see e.g. [22]), too high to solve the strong CP problem. Model-dependent 
axions arise from the internal components of the B— field. Unfortunately, their decay con- 
stants are in general not much lower than in the previous case [22]. Admissible axions, 
on the other hand, might arise from the breaking of accidental global U(l) symmetries 
realized as low energy remnants of (stringy) discrete symmetries. In this case, the axion 
decay constant shall depend on the VEVs of the fields responsible for the breakdown of 
such global U(l)s. 

In general, there is a set of fields with a large VEV dictated by the FI term (one or 
two orders of magnitude below the string scale). This is the scale where (some of) the 
gauge U(l) symmetries are broken and where also some of the anomalous U(l)s will be 
broken. The latter ones should not have a QCD anomaly (they could have e.g. hidden 
sector anomalies) and/or should be candidates for a quintessential axion. There could in 
addition be another set of singlet fields that obtain smaller VEVs that break the relevant 
PQ-symmetry at a scale in the axion window and thus give a hierarchically smaller F a . 
The requirements for a satisfactory model are thus clear: 

• find a model with an accidental (color)-anomalous U(l)*, 

• identify a vacuum configuration where the VEVs driven by the FI term do not 
break U(l)*, 

• search for a vacuum configuration where U(l)* is broken by a VEV in the axion 
window (some other gauge U(l)s might be broken here as well), 

• check that higher order non-renormalizable terms that break U(l)* explicitly are 
sufficiently suppressed to avoid a too "large" axion mass. 
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A comment regarding the latter requirement is in order. Higher order terms in the 
superpotential break explicitly the accidental U(l)* symmetry. This can pull the axion 
off from the CP-conserving vacuum and generate an additional contribution to the axion 
mass m\ [23,24]. m\ is constrained by the bounds on the electric dipole moment of the 
neutron [25,26]: 

ml < 10- 11 m% CD « ID" 15 ^j! , (17) 

"a 

where m 2 Q CD 10 -4 GeV 4 /F% is the axion mass induced by QCD instantons. Not 
every U(l)*-violating superpotential term, however, contributes to the axion mass. The 
magnitude of the axion mass, induced by the explicit breaking of the U(l)* symmetry, 
does not necessarily depend on the order of the lowest term in the superpotential that 
breaks this symmetry. It rather depends on the order of the (lowest) U(l)*— violating 
term that contributes to the axion mass directly. Note e.g. that if a superpotential term 
breaks U(l)* explicitly, but contains two or more fields that do not develop VEVs, it 
does not contribute to the axion mass. It follows then that verifying whether eq. (|17|) 
is satisfied depends strongly on the specifics of the vacuum configuration and must be 
analyzed case by case. 

We shall explore some candidate models in the next section. 



3 QCD Axions in Orbifold Models 

Let us turn now to the question of whether we can find suitable QCD axions in heterotic 
orbifold models. Recent progress in this framework has revealed that models resembling 
many features of the MSSM can be constructed [27-31]. In particular, we analyze a 
subset of the models found in [15] and identify those with the following properties: a) 
their superpotentials W must be symmetric under an accidental PQ symmetry, and b) 
they should admit vacua such that the PQ symmetry gets spontaneously broken at a 
phenomenologically admissible scale (see eq. (P)). 

Identifying PQ symmetries of the superpotential W of a model requires to verify 
whether all superpotential terms of a given order preserve some U(l) symmetry other 
than the gauged ones. To find these accidental symmetries, it is convenient to use the 
method of [14], which we briefly summarize in the following. The superpotential of a 
model can be written as 

# terms # fields 

W= W with Wi= Yl (18) 

» 3 

where all coefficients are set to one and A^- are non-negative integers constrained by the 
order n of the superpotential according to Y2j ^ij — n i f° r an Each element of the 
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kernel (or null space) of the matrix (A T A) defines the charges of the fields under an 
Abelian symmetry of W. Subtracting those symmetries contained in the gauge group, 
one is left with all accidental symmetries denoted by U(l) acc . 

Not every accidental symmetry of the theory can play the role of a U(1) P q. It is 
crucial that such symmetries also produce mixed U(l) acc — SU(3) — SU(3) anomalies, 
that is, 

^ = X>^cc^0, (19) 

3 

where qi cc denote the accidental charges of the fields (fj and lj refer, as before, to the 
quadratic indices of the SU(3) representations. 

Using the method outlined above, we search for the accidental symmetries in a subset 
of realistic Zg-II orbifold models found previously in the so-called Mini-Landscape [15, 
17,18]. We use the notation introduced there. We analyze 55 orbifold models based on 
the Mini-Landscape shift ySCK 10 )- 1 exhibiting the following properties: 

• SM gauge group, 

• non-anomalous hypercharge U(l)y C SU(5), 

• spectrum = 3 generations + vector-like exotics, 

• heavy top, and 

• super symmetric vacua. 

The last property deserves a comment. The pseudo- anomalous symmetry U(1)a 
(generically present in orbifold models) induces the FI term [32] 

M 2 

Da « E^K^+ig^lEd- (20) 

3 3 

Therefore, in supersymmetric vacua a set of fields ipj have to acquire non-vanishing 
VEVs, such that Dj± = 0. In general, solutions of the D = conditions can be expressed 
by gauge invariant monomials [33], 

/ = <vr---> (2i) 

where mi are positive integers. In the case of an anomalous U(l), such a monomial 
must be invariant only with respect to the non-anomalous gauge symmetries in order 
to cancel the FI term, eq. (pOj) . Since in our conventions ^ ■ q J A > 0, a suitable vacuum 
configuration exists if there is a set of SM singlets, Si C ipj, that build a monomial with 
a total negative anomalous charge, ^ ■ rrij qj^ < [30]. Therefore, the condition D\ = 
sets the VEVs of the singlets Si and, thereby, the scale of the breakdown of U(1)a- In 
heterotic orbifolds, this scale can be as high as ~ 10 17 GeV. All models studied here 
have vacua satisfying D = and preserving the SM gauge group unbroken (including 
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Order of W 


# Models with some 

U(l)acc 


# Models with some 

U(1) PQ 


3 


55 


55 


4 


55 


54 


5 


55 


15 


6 


34 






Table 1: Statistics of frequency of accidental PQ symmetries in a subset of promis- 
ing ^6-11 heterotic orbifold models. 

hypercharge). [§ This contradicts recent statements about model building in heterotic 
string orbifolds [34]. 

For each model, we choose a monomial suitable to cancel the FI term that involves 
only a few SM singlets. The VEVs of the fields in such a monomial is set close to the string 
scale by eq. (120p . Further VEVs for other SM singlets which we will consider later are 
assumed to be significantly smaller. The global U(l)s in which we are interested are those 
of the effective superpotential at finite orders obtained after the cancelation of the FI 
term. We find that the studied orbifold models have generically one or more accidental 
symmetries. As presented in table [IJ although most of the models have approximate 
global symmetries up to high orders, only at orders 3 and 4 most of them exhibit U(l) acc — 
SU(3)— SU(3) anomalies. In fact, at order 6 in W, no model of the limited sample analyzed 
in this study has symmetries which can serve as a U(1) P q. 

The explicit breaking of all accidental PQ symmetries at order 6 in these models 
could spoil the CP-conserving vacuum. However, from the discussion of sec. 12.41 a solu- 
tion to the strong CP problem is still possible provided that the PQ-violating couplings 
contributing to the axion mass appear at higher orders. In addition, we have seen that 
requiring a hierarchy of VEVs lets us set the axion decay constant at an intermediate 
scale. Such a hierarchy could arise in this type of models from other symmetries of the 
low energy theory [13]. In the following we examine a model with only one U(1) P q. 



3.1 A Z 6 -II Model with a QCD Axion 

A 2^-11 orbifold model which can lead to a suitable QCD axion is defined by the following 
gauge shift and Wilson lines: 

v«xu)a = fi-i-ioooooV-- 1 --- 1 --- 1 --- 

v 3' 2' 2 ' ' ' ' '7 V 2 6 ' 2 ' 2 ' 2 ' 2' 2' 2 
/ 1 111111 1\/ 1 22 5 21 

Wz ~ l _ 2' _ 2'6' 6' 6' 6' 6' 6/ I 1 ' ' - 3 ,_ 3' 3'~3'~3'3 



3 As we only know the form of the superpotential without coefficients, we have to assume that F — 
can be fulfilled and stabilizes the remaining D-ftak directions. For generic coefficients of order f, this has 
been checked in some examples. 
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# 


Representation 


Label 


# 


(Anti-)Rcpr. Label 


3 


(3,2 


M)l/6 


Ii 






9 


(1,2 


l,l)-l/2 


ii 


6 


(l,2;l,l) 1/2 ii 


3 


(1,1 


l,l)l 


ei 






3 


(3,1 


l,l)-2/3 


Ui 






8 


(3,1 


M)l/3 


di 


5 


(3,l;l,l)_ 1/3 di 


12 


(1,1 


l,l)l/2 


4 


12 


(l,l;l,l)_ 1/3 sr 


2 


(1,1 


l,2)i /2 


4 


2 


(l,l;l,2)_ 1/2 xr 



# 


Representation 


Label 


4 


(1,2;1,1) 


mi 


2 


(1,2;1,2) 


Ui 


65 


(i,i;i,i)o 




12 


(i,i;i,2) 


hi 


10 


(1,1;4,1) 


Wi 


10 


(i,i;4,i) 


Wi 



Table 2: Massless spectrum. The quantum numbers are shown with respect to 
SU(3)c x SU(2) L x SU(4) x SU(2), the hypercharge is given by the subscript. 



W 2 



3 1111 1 1 1 



5 3 _5 _5 

2' 2' ' 2' 2' ' ' 



The unbroken gauge group after compactification is 

g = SU(3)c7 x SU(2) L x U(l)y x [SU(4) x SU(2) x U(1) A x U(l) 7 ] 



(22) 



(23) 



where U(l)y denotes the standard SU(5) hypercharge and U(1)a, the anomalous U(l). 
The massless matter spectrum is given in table it contains three MSSM generations 
plus vector-like exotics with respect to SU(3)c X SU(2)l x U(l)y. All non-Abelian and 
Abelian charges of the massless particles are provided at our web page [35]. 

The selected vacuum has to preserve N = 1 SUSY. With this purpose, we require 
that the fields contained in the gauge invariant monomial 



Ifi 



,0 \ 2 _0 / \3 / ^ 2 



(24) 



(with total U(1)a charge <?a(^fi) = —25/6) acquire VEVs such that the FI term is 
canceled. As a consequence, the gauge symmetry is further broken to 



Q' = SU(3) C x SU(2) L x U(l)y x [SU(4) x SU(2) x U(1) E 



(25) 



After symmetry breaking, the effective superpotential W of this model contains 870 
couplings of order up to five and is invariant with respect to 19 accidental U(l)s, in 
addition to the gauge symmetry group g' . Only one accidental U(l), that we denote 
U(1) P q, presents mixed PQ-SU(3)c anomalies: 



If the fields entering the gauge invariant monomial 



(26) 



T _ „0 
1 — s l s 36 s 37 



(27) 
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Til 

Label 


9A 


1Y 


93 


94 


95 


96 


97 


98 


99 


9PQ 


„o 
b i 


5 
3 





-1 


-3 








-10 








-2 


„o 


11 





i 


ii 


5 


1 


1 


i 


2 





°19 


18 


2 


6 


3 


2 


2 


3 


3 


„0 
6 36 


22 
9 








4 








2 











„0 
A 37 


13 
9 





1 


-1 








8 








1 


S° 
6 49 


7 
9 





1 

2 


4 
3 


5 
6 





12 


2 
3 


4 
3 





„0 
6 52 


22 
9 








-4 








-2 











„0 
6 57 


4 
9 





1 

2 


1 

3 


5 
6 





-14 


2 
3 


i 

3 





„0 


13 





1 


10 


-> 


1 


9 


1 


2 





*59 


9 


2 


3 


6 


2 


2 


3 


3 



Table 3: U(l) charges of the relevant non-Abelian singlets of a Z^-II orbifold model 
with a PQ symmetry. 

acquire VEVs, U(1) P q is broken spontaneously and the axion 

a = M ( 9 PQ Viai + 9PQ U36a36 + 9PQ U 37037) (28) 

appears in the theory. Here aj, qpQ and Vi denote respectively the phase components, 
PQ charges and VEVs of the singlets of eq. (|27p . The U(1) P q breakdown scale M is read 
off from eq. ©. 

Along the D-flat direction described by eq. (|27p . the VEVs of the singlets are equal. 
Thus, the axion decay constant is given by 




where we have considered the singlet charges displayed in table O Provided that the 
VEVs Vi are well below the scale of the FI term, say at ~ 10 12 GeV, the inverse axion 
coupling can agree with the current experimental bounds, eq. ([1]). Note that we have 
disregarded the axion-like components of the fields that acquire VEVs at the larger scale 
(see eq. (|24p ) as their contributions are subdominant. 

There are two issues that still have to be addressed. First, since the U(1) P q is broken 
at order 6, we have to verify whether there appear large contributions to the axion mass. 
In the current model, the first contribution to the axion mass from the explicit breaking 
of U(1) P q arises at order 16. We have estimated the contribution to the axion mass 
to be m 2 ~ 10 -18 GeV 4 /.F 2 , consistent with the current experimental bound, eq. (fT7|) . 
Second, in the vacuum configuration chosen here, most of the exotics remain massless. 
This problem is due to the small number of fields that are chosen to develop large 
VEVs (see eq. (|24p ) . However, as shown in [15-19], there exist supersymmetric vacua 
where all exotics acquire masses because (almost) all the SM singlets attain VEVs. In 
the current model, we find a monomial Ifi containing 32 fields that leaves a U(1) P q 
symmetry unbroken and provides masses for many more exotics. Whether there exist 
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configurations where all exotics are decoupled and simultaneously a U(1) P q is retained 
is a question beyond the scope of the current letter and will be addressed elsewhere. Note 
that this is related to the magnitude of m\ , as integrating out heavy exotics could induce 
further contributions to the axion mass. 

4 Conclusions 

Discrete symmetries (as found abundantly in models of the heterotic brane world) are 
suited to solve the strong CP problem in a satisfactory way. They might give rise to 
accidental global U(1) P q symmetries as "low energy" remnants. Accions - the axion 
like particles of the spontaneous breakdown of these symmetries - could have a decay 
constant F a in the allowed axion window, eq. (H|). 

The key observation for this latter fact has been explained in section [21 formula (]14p . 
In a system with an anomalous and a non-anomalous U(l) symmetry broken by two fields 
with nontrivial VEVs, it is the smallest of these VEVs that corresponds to the accion 
decay constant F a . This allows a decoupling of F a from the large VEVs required to 
cancel the FI terms (which appear to be only one or two orders of magnitude below the 
string scale). The actual value of F a depends strongly on the vacuum configuration of the 
system and requires further model building. In this framework of models with accidental 
U(l) symmetries, such a hierarchy of the scales could be potentially created along the 
lines explained in ref. [13]. 

The implementation of the accion-system in the framework of realistic MSSM models 
from the heterotic string is still in its infancy. The identification of suitable candidates 
for accidental U(l)'s requires an enormous amount of computer work. In section [3] we 
have examined a few candidate models from our previous MSSM search that can serve as 
an existence proof for the accion solution to the strong CP-problem. A few things remain 
to be done: the decoupling of exotics at higher order in the superpotential and a detailed 
discussion of the accion mass arising from the explicit breakdown of the accidental U(l) 
symmetry. Such an analysis is beyond the scope of the present paper. It would require 
substantially higher computer power and/or new algorithms to efficiently attack the 
problem. Work along this direction is under way. 
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